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Abstract. White-light images are presently the primary source of information on phvsical conditions 
m the solar corona at distances greater than a few tenths of a solar radius above the limb. As a 
consequence, we still only have an incomplete description of structures extending beyond the solar 
hmb. In particular, streamers, although observed for decades, represent a poorly known phenomenon. 
SOHO, to be launched in 1995. will be able to make long-term observations of these features up to 
heights ot a tew tf . . both in white light and UV. In this paper we present simulations of La intensitv 
iooi r0na strcarners " b ase d on the two-dimensional t2-D) model developed by Wane et al ( 1992 
1993) via a time-dependent numerical relaxation approach. Because the model is 2-D we make an 
a priori hypothesis about the extension of streamers in the th.rd dimension. La data, obtained from 
a rocket (Kohl et. al.. 1983). allowed us to identify a shape which tits the observations. We consider 
streamers with different magnetic held contigurations and at different position angles with respect 
to the plane ot the sky to illustrate how different regions along the line of sight contribute to the 
emergent intensity. Our purpose is twofold: to provide guidelines for UVCS observational operations 
and to explore the parameter space in order to understand the role of geometric factors and of the 
physical state ot the corona in determining the overall streamer brightness. We conclude by showing 
how the results guide the tuture development of streamer models. 


1. Introduction 


Before the advent of coronagraphs, eclipses offered the only means to observe 
the solar corona. In spite of the short time over which coronal structures were 
visible, their basic characteristics have been reproduced in beautiful drawings (see. 
e g., Foukal, 1990), which prove that streamers - the most prominent white-light 
coronal features - had been known and observed for decades. Following earlier 
balloon-borne coronagraph experiments (Newkirk and Bohlin, 1965; Dollfus, Fort, 
and Morel, 1 968), different techniques have become available, and radio telescopes 
(see, e.g., Gopalswamy, Kundu, and Szabo, 1 987) and space experiments (see, e g 
for Skylab, Hildner et al, 1975; Poland, 1978; for SMM, Illing and Hundhausen, 
1 986; Kahler, 1991) collected a wealth of data on these large-scale features, which,' 
nevertheless, remain poorly known. We know that streamers, at the time of solar 
activity minimum, are concentrated along the solar equator in a belt, which broadens 
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During past several months we have established a two-dimensinal time-dependent 
numerical MHD model to do the numerical simulation of large scale coronal structure with 
coronal streamers and coronal holes. The motivation for this study is to develop a tool for 
the interpretation of the observations to be obtained from UVCS and LASCO experiments 
on board the upcoming SOHO mission. 

First, we have simulated coronal streamers and given four examples. Please see at- 
tached papers for the details. 

Then we simulated the streamers and holes simultaneously. In this study the ini- 
tial density and temperature given are not homogeneous from the pole to the equator. 
The density decreases from the equator to the pole. The temperature increases from the 
equator to the pole to get high velocity at the pole, since we used the polytropic state 
equation instead of the real energy equation and this temperature is effective temperature. 
According to density and temperature distribution the initial Parker’s solution has been 
obtained along each altitude direction from the pole to the equator. In the steady state at 
the hole region the high solar wind with low density has be achieved, and at the streamer 
region the low solar wind with high density, especily at the closed field region with almost 
zero velocity and high density, has been achieved simultaneously. The results are entitled 
’’Simulation of Streamers with Coronal Holes” by S. T. Suess, S. T. Wu, A. H. Wang 
and G. Poletto were presented at the Second SOHO Workshop at Elba, Italy and will be 
published in the proceedings (in press) Kluwer Publishers, 1994 . 
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to cover a wider range of latitudes, at other epochs. We know that streamers are 
stable structures, which may last for several rotations (see. e.g., Poland, 1978). 
However, we know little about the formation, or the disruption, of streamers: we 
have only a few observations, for instance, of the contraction of a post/mass ejection 
structure to form a streamer (see, e.g., Illing and Hundhausen, 1986). As to the 
physical properties of streamers, our knowledge is limited, as we will discuss later 
on, to an incomplete understanding of their density structure. This situation may in 
part be ascribed to the tact that spectroscopic instalments were unable to operate 
at the large heights reached by streamers, which therefore have been observed 
mainly in white light. On the other hand, the solar community seemed to have little 
interest in streamers: no review article has been dedicated to these structures for 
years, notwithstanding the steady proliferation of scientific papers. 

However, streamers are relevant to a number of problems. They represent 
the ideal structures to investigate the differences in temperature, density, flow 
velocities, and magnetic held structure between magnetically closed and open areas 
and to gam some insight into the physical conditions of current sheet regions, which 
purportedly occur in streamers. Koutchmy (1988) pointing out how tangential 
discontinuities are almost systematically parts of large streamers, noticed how 
their analysis would offer the best estimate of the magnetic held value in the 
corona, provided that the temperature keeps constant across the discontinuity. The 
capabilities of streamers in providing data crucial to the solution of these questions 
have hardly been exploited. 

In solar wind physics streamers plav an ambiguous role. It has been claimed (see, 
e.g., Feldman et al ., 1981) that they constitute the major source of the interstream 
and low speed solar wind, but we know neither how large is the contribution 
of streamers to the solar wind mass llux, nor the geometry of the open field 
lines associated with streamers, along which the wind purportedly propagates. The 
reason for the depletion of helium, which seems to accompany the slow wind from 
streamers (see, e.g.. Gosling et al., 1981), is not well understood. 

In the future, SOHO instrumentation will offer us a means to learn more about 
these structures. In particular, UVCS will be able to make EUV observations of 
streamers, up to < 1 0 over an extended period of time, thus allowing us to get a 
new kind of data whose capabilities have not yet been explored. In order to provide 
guidelines to be used in devising UVCS observational sequences, we present, in 
this paper, a variety of simulated, typical La observations of streamers. As a basis 
for our simulations, we adopt the two-dimensional streamer model, developed by 
Wang et al, (1992, 1993), which is summarized in Section 2. After comparing, in 
Section 3, the model predictions with observations of densities in streamers, we 
calculate, in Section 4, the La emission from different streamer configurations and, 
in Section 5, we simulate La observations of a streamer carried around by solar 
rotation. Finally, in the Discussion, we illustrate some future development of our 
work. 
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Of the magnetohydrodynamic equations for two-d,mens,ona .' Z s^r^nde 

The hn h simulation is further described by Wan 2 et al. ( 1 999 , 991) 

The boundary conditions at 1 R,. are that the temperature and density are con 
s ant m latitude and that the vector magnetic held is a potential magnetic held 
Three magnetic held geometries are used: a dipole, a qualpole. aLTh^o 
the scalar potentials are proportional to Afcosd). A(costf) and PAc os^f r ’ 
spectively. There are two dimensionless numbers: the polymop,^ £ 



* W UI1U .y — : ' 'S 

temperature and density are 1 .8 x 1 0 6 K and 25 > 
they are 1.44 x l0 6 Kanri5 6l v I0 7 cm’ 3 - ■ 


1 0 8 cm 3 . For the low ti case. 
The three magnetic held geometries 


they are 1.44 x I0 6 K and 5 61 x 
naturally lead to a single equator, alM^er 

respect?vdy “ m ‘ d ' lat,tUde Streamer for the dl P° le - q^dmpole. and hexapole. 

honnHn Simul “ U ? n extends from 10 77 to > 5.0 R, and from pole to equator The 
boundary conditions are symmetric about the equator, so a solution in the opposite 

hemisphere is not necessary. In this paper, we only quote results inside 7 0 R 
because th,s covers the range over which closed streamer stmctums mo fLuen^v 

Tr; 7 0 2 °„ 8rid r mst ™ ,hei -«^“S3S£ 

■ , .. nd 7.0 R^. The initial state consists of a potential held and the 

so ution or a spherical.y-symmetric wmd for the given base'tempemmm dens tv 

““ — • *-*• a"d velLy JZ* 

nown in figure I. The temperature curves appear irregular due to the small 
change temperature over the relatively large radial range - a consequence of the 

£ ca7utatTo„ e so be 'h n8 , " ear ° nly ' hree S ' gniNcam ^ "Trained after 
the calculation so what is seen here is essentially roundoff error. 

Kesuits for the steady-state solutions, given the above initial conditions are 
ottered according to the four cases treated: (a) dipole, i = |. 0 , !£%££ 

field L°’ } heXapole ’ = 1 and (d > di P° le ' J = 0.2. The steady-state magnetic 
held geometries are shown in Figure 2. Here is seen the well-known property that 

ow is nearly radial beyond 3-4 The streamers are those volumes which are 
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Radius (solar radii) 



Radius (solar radii) 


Fig. la-b. Density, temperature, and velocity protiles in radius of the initial (£ = 0) state in the 
relaxation, (a) Profiles for the d = 1.0 cases, (b) Protiles lor the d = 0.2 case. Note that, except 
for the velocity scales, the scales differ between the two panels. The irregularities in the temperature 
protiles are roundoff errors. 
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magnetically closed, and it is evident that relatively small volumes in the streamers 
remain magnetically closed in comparison to the initial states. These volumes are 
surrounded by a low-density shell but. as will be shown below, the densities in the 
arge coronal hole-like open regions are otherwise only slightly lower than in the 
streamers. In each panel of Figure 2, four dotted lines are shown and labelled A 

B, C, or D These lines indicate the positions used to plot variables versus radius 
in Figures 3 and 4 below. 

In displaying results for the physical variables, we will concentrate on the density 
and the velocity, these being the two variables that determine the La intensity. Other 
variables will only be shown when necessary for physical understanding. Figure 3 
shows the density variation in the radial directions labelled in Fieure 2 These 
plots show the relative density - the density divided by the density pFotted in 
Figure I . Thus, it measures the change in density with respect to a spherical How 
profile. Curves D in Figures 3fa). 3(c). and 3(d), and curve C in Figure 3(b) are 
all profiles cutting through the cores of streamers. It is seen that there is a density 
enhancement in the core ranging from =75% for the low d dipole down to ?0% 
tor the equatorial streamer in the hexapole. On the Hanks of streamers, for example 
as shown by curve C in Figure 1(a), there is a density deficit. Nevertheless in the 
centers ot open regions, the density deficit is always less than =20%. Thi's is an 
important point to note in applying this specific model to computing La intensities, 
and we will return to it later. Also apparent in this Hgure, and the ones that follow 
s owing variation of other variables with radius, is a rapid fluctuation from one 
grid point to the next for the first two or three points above the base. This is a 
consequence ot the type of extrapolation used to determine those variables at the 
boundary which are calculated from the interior solution instead of fixed by the 
boundary conditions (Wang et «/., 1992, 1993; Steinolfson, Suess, and Wu. 1982). 
The phenomenon has no sigmHcant effect on the solution above the fourth grid 

point and, in particular, no effect above 1 .2 R.„ - the minimum radius UVCS can 
observe. 

Figure 4 shows the radial velocity along the directions labelled in Figure 2. In 
addition, the initial state prorile is plotted as a dashed line - the same profile as 
s own in Figure I. The How speed is. except for numerical diffusion, identically 
zero in the streamers. This is obvious from curves D in Figures 4(a), 4(c), and 4(d), 
and curve C in Figure 4(b). Above the streamers, the velocity is greatly reduced 
However, away from the streamers, whether in the center or at the edges of open 
regions, the How speed differs little from the initial state. 

Taken together, the density and How speed illustrate that this is not a constant 
mass Hux model - the flow speed at the base varies with polar angle. Therefore, 
these results appear somewhat different than might be anticipated based on intuition 
gained from, for example, the calculation by Kopp and Holzer (1976). A rapidly 
diverging magnetic field does not necessarily lead to a low density and high velocity 
flow - depending also on whether the rapid divergence occurs below or above the 
some critical point. So. on the Hanks of streamers, where rapid field line divergence 


</. N( )C I ET AL 


(b) Quadrupole 


liexapole 


Dipole. 0 = 0 2 


Fig. 2a-d. Magnetic field line plots for the four cases: (a) dipole. 3 = 1.0, (b) quadrupole, (3 = 1.0, 
(c) hexapole, 3 = I 0, and (d) dipole. 3 = 0.2. The relaxation times allowed to reach these equilibria 
are. (a) 22.22 hours, (b) 16.67 hours, (c) IS. 06 hours, (d) 19.44 hours. In each panel, four dotted 
lines arc labelled ‘A. B. C. or D\ These show the radial directions used for plotting certain variables 
vs radius in succeeding figures. Thus, the quadrupole in plot <b) will have these variables plotted vs 
radius at the pole (A), at the edge of the polar open region i B). through the mid-latitude streamer (C). 
and along the middle of the equatorial open reeion ( D). 
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?h g ,' DCn ^‘ y " S a ‘ UnCti0n °‘ r3dius ' Each P anel IS for lhe corresponding case in Figure 2 

Soiled A “b C°n ’ i ‘ ,C d h lre t C1, ° n Y hOWn m RgUrc 2 - For exam P le - four curves in (c). 
n • a [ C 3 S thC l0Ur dlrccIlons sh °*" ■" Figure 2fc) and labelled in the same 

rZ? T m 71 r UrVC h3Ve bCen normalizcd [0 the,r 'he initial profile (Figure 1(a) 

thin nd Figure 1(b) tor case (dll. Hence, a density enhancement is indicated by values greater 
uni v, and vice versa. The density concentrations in the streamers are clearly visible generally 
bemg on the order of 25* to 50% above the initial state. generally 



occurs over a narrow range, the density does become small (as in Figure 3(a), 
curve C) but the flow speed never gets large (as in Figure 4(a), curve C). The 
other important thing to note is that the cusp (top of the streamer) in the j3 = 0.2 
dipole lies at about 6 This seems large enough to cover the range of streamer 
heights expected in the solar corona and therefore we will not concern ourselves 
with computing models for smaller d-values. 

A more complete picture of the behavior of the density can be gained by also 
considering plots of the density versus polar angle at different heights. These are 
shown in Figure 5, where each curve is labelled with the heliocentric distance it 
represents (e.g., 2.30 R, is at 2.30 R . ). Figure 5(a), the t) = 1 .0 dipole, shows the 
density enhancement in the streamer (polar angle of 90°), the deficit in the adjacent 
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(c) 


Hcxopole 



(d) 


Dipole, t3 = 0.2 



Fig. 4a-d Radial velocity as a function of radius. Each panel is for the corresponding case in 

,'v' e P '° tled al0ng the direclions shown ln Figure 2. For example, the four curves 

‘ abelled 1 A. B. C. D . arc along the lour directions shown in Figure 2(c) and labelled in the 
same manner. The dotted line in these plots is the t = 0 profile, the same as shown in Figures 1(a) 


trough, and the large plateau of density that is only a small amount less than in 
the streamer and extending throughout the open region beyond the trough. The 
behavior of the density around the mid-latitude streamer in the quadruple, and 
around the mid-latitude and equatorial streamers in the hexapole is very similar 

The only difference for the 6 = 0.2 dipole is that the troughs are considerably 
broader. 

The broad, high density plateau in the open region is distinctly unlike a coronal 
hole. The reason for this is that in this model no effort has been made to generate 
the high coronal hole flow speeds that lead to low densities. Suess et al. (1977) 
have shown in a similar model that a temperature increase of 50% or more at the 
center of the open region is necessary in a polytropic model such as this to produce 
densities like those that are observed in coronal holes. 
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( a ) Dipole 



(b) Quadrupole 



Polar angle (degrees) 


De " su y vs polar angle, between the pole <0") and the equator (90°). Each of the curves is 
“belled according „ lhe hcl,ocemn c distance it refers to. Thus, the curve labelled I 70 R indicates 
£! ? 1 •» J* *• °““™ The d.„ s ,„ „ b*. i s cons™ and „ 

tronoh in ^ . hC b ^ se> ,here . ls a sma “ densil y enhancement in the streamer (ca. 5% to 50%) and a 

rough in density at the edge ol the streamer. In the middle of the open region, the density is very close 

° “ WaS m "'‘"h St3te ,SCe alS ° FigUreS **** « is no. an orde^fTagnnudl 

the hnie t ! an observed coronal hole, is that we have used constant temperature and density at 
reauired \n° Pr ° dUCe 3 coronal ' hole - |lke P ro,ile a polytropic model such as this would have 
oSw.se ^tated7= 1 .0 l ' mpemUK 31 the base °' the ^ re g'°" (Suess et al.. 1977). Unless 


Figure 6 shows the radial velocity plotted in the same manner as the density in 
igure 5. The velocity is again seen to be essentially zero inside the streamers), 
w ose height decreases rapidly with increasing magnetic field complexity Thus 
while the 8 = 1 .0 dipole streamer extends to 3 neither of the 0=1.0 hexapole 

streamers reaches beyond 1 .70 R,-,. As indicated above, the flow speed throughout 
the open region is very similar to the initial state flow speed, excepttng for small 
humps on the flanks of the streamers. 

These four models constitute the basis for the calculation of La intensities. 
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^ Dipole 



(k) Quadmpolc 



Polar angle (degrees) 


rofnSwn, h a r al vcloc, ' y versus P° lar an glc. between the pole and the equator. Each curve, plotted 
them ? m locentnc distance. is labelled in the same manner as in Figures 5f a— d) The velocity in 
the magnetically closed regions is essentially zero. The reason it is not identically zero ts that there is 
a smal amount ol t numerical diffusion - quite small, as indicated by the velocity^ess than 10 km s' 1 
ms.de the 3 = I 0 dipole streamer at 2.30 A, (top left panel). 


In using them, we progress from a stra.ghtforward calculation of the intensity 
measured when v.ewmg the streamer as seen in Figure 2 (i.e., from a pos.tion in 
t e magnetic equatorial plane) to other viewing positions and to approximations 
based on the models. A comparison with intensities from observed streamers 
is used to guide suggestions for further development of the model by illustrating 
specific weaknesses in the present four models. We will conclude that a satisfactory 
physical model of streamers, for the purpose of computing expected UVCS La 
intensities, can be constructed through the application of the present simulation with 
an appropriate choice of boundary conditions to better represent the dynamics of 
the solar wind in the open magnetic held regions. This is well within the capabilities 
ot the simulation and will constitute the next stage of this project. 
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As mentioned in the Introduction, density is the only physical quantity in stream- 
ers to be even partially measured. Coronagraph images usually are polarization 
brightness i pB) images. Because pB is proportional to the line-of-sight intearal of 
the density times a scattering function (Billings, 1966), it is possible from pB im- 
ages to derive density maps. Usually, this procedure is performed with the Van de 

Hulst method ( 1 950): that is. the density distribution is supposed to be cylindrically 
symmetric. 

This technique has been used to derive densities from eclipse observations 
of streamers. Dollfus. Laffineur. and Mouradian (1974) derived densities for a 
number of streamers observed in the eclipse of February 15, 1961 . and compared 
their values with those pertaining to 13 different streamers, which represent all 
results published between 1 952 and 1 972. From this, as well as from a comparison 
between values derived by different authors for the same streamer, these authors 
conclude that different determinations for the same structure agree only within 
a factor 2 and different structures may have densities which differ, at the same 
altitude, by a factor 10. even in streamers observed at the same eclipse, that is, 
independently of the epoch of the solar cycle. All densities refer to the streamer 
axis, Dollfus < t nl. assumed that streamers are axially symmetric and that the 
distribution of density, in the direction normal to the streamer's axis, is somehow 
intermediate between being uniform and having a gaussian distribution. Different 
assumptions on the streamer geometry, or on the distribution of density across a 
streamer, may possibly explain some of the discrepancies in the values derived for 
the same structure. 

Densities predicted by a theoretical model have to comply with this rather loose 
observational constraint. Figure 7 gives, on the left panel, the behavior of the 
density predicted by our model along the axis of equatorial streamers, for dipolar 
and hexapolar configurations in the case of .) = 1 .0, and for dipolar geometry only 
in the low 6 = 0.2) case. Densities along the axis of off-equator streamers, both 

in the quadrupoiar and hexapolar geometries, are approximately equal to those 
along the axis of equatorial streamers and are not shown. In the right panel, we 
present a figure made up from Figures 16(a-b) of Dollfus et al . , which shows, 
besides all density determinations in streamers between 1952 and 1972, the values 

derived by Dollfus et al. for four different streamers observed in the eclipse of 
1961. 

The J = 1.0 curves, cutting through the bundles of curves shown in the right 
panel, represent correctly the observed densities. As we said, it is not possible to 
establish, from the data published so far, any trend in the different behavior of 
individual streamers. For instance, contrary to expectations, the streamer closest 
to the equator in Dollfus et al. s data (position angle 95°), has lower densities 
than a mid-latitude streamer (position angle 55 3 ) and the highest densities among 
those from published results 1952-1972 pertain to a streamer observed close to 
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Radius (solar radii) 

Fig. 7. Left: density vs height proliles along the axis of equatorial streamers in dipolar ( 3 = 1.0 
- solid line - and ii = 0.2) and hexapolar geometries as derived from Wang et al.'s ( 1992. 1993) 
models. Right: densities along streamers from all results published between 1952 and 1972 (labels 
1-13. solid lines) and along 4 streamers observed at the February 15. 1961 eclipse (dashed lines, 
labelled by their position angles). See Figures I6(a— b) ot Dolltus. Latfineur. and Mouradian ( 1974) 
for further information on this panel. 

the activity minimum (curve 4, February 5, 1962 eclipse). Hence, we can only 
conclude that the 3 = 1.0 curves are consistent with observed density profiles 
in coronal streamers. This conclusion is confirmed by a comparison of the radial 
density distribution predicted by our model with densities derived from Clark Lake 
Radioheliograph streamer observations (Gopalswamy, Kundu, and Szabo, 1987). 
This comparison is limited to the lower corona, at heights below and near 2 Rq, and 
shows (Figure 8) how our predicted density profile lies between densities derived 
from fundamental and harmonic plasma hypotheses. Finally, we notice that the 
lack of a definite observational difference between low and high latitude streamers, 
agrees with predictions from our model. 

On the other hand, by comparing the two panels of Figure 7, we conclude that 
our 3 = 0.2 case is not realistic because densities are too low and the cusp is far 
too high. Nevertheless, the 3 = 0.2 curve show's a marked change in its slope that 
is not so evident in the 3 = 1 .0 hexapote streamer, and is altogether absent in 
the 3 = 1.0 dipole streamer, which may reproduce the behavior shown by some 
of the observed structures. To recover this break in the density profile at higher 3 
probably requires changing conditions outside the streamer. A change in the density 
gradient of the observed profiles has, in fact, been interpreted in terms of a different 
behavior of this physical parameter in the region of the streamer’s helmets (Dollfus 
et at., 1974). Obviously, the present simulations do notallow us to predict whether 
the resulting curve will be capable of reproducing some of the observed density 
profiles more closely than the high 3 curve. However, it is likely that structures 
with differing cusp heights correspond to different 3 values (Steinolfson, Suess, 
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Fig. 8. Predicted density vs height prohle alone the axis of the cquatonal streamer in the case of hieh 
■ i (i) = 1.0) dipolar configuration (solid line) and densitv determination from radio observations of 
a streamer obtained on July 27. September 12 and 17. 1985: o represent values of density derived 
Irom the hypothesis ot fundamental plasma emission. A represent values of densitv derived from the 
hypothesis of harmonic plasma emission. 


and Wu, 1982). We conclude that the comparison between model-predicted and 
observed density profiles points to the need for a thorough analysis of the effect 

of different boundary conditions in and around streamers on the resulting densitv 
profiles. 7 

In the following section we proceed to evaluate the La emission in streamers 
adopting the 3 = 1 .0 models. 


4. La Emission from Streamers 

The formation of the La line in the solar corona has been discussed by a number of 
authors (Gabriel, 1971 ; Beckers and Chipman, 1974; Withbroe et al., 1982; Noci, 
Kohl, and Withbroe, 1987) who showed how coronal La observations can be 
used as a diagnostic tool to determine coronal densities, temperatures, and outflow 
velocities. Although, at coronal temperatures, only =sl proton in 10 7 is tied up in 
neutral hydrogen, the strongest component of the coronal La is due to the scattering 
ot chromospheric La photons by neutral hydrogen atoms. An electron scattered 
component, produced by Thomson scattering of La radiation, is about three orders 
of magnitude weaker than the resonantly scattered component and will be ignored 
in the following. 

The total (i.e., integrated over the line profile) La intensity, as observed along 
the direction n is given by 


/ = 


h B } 2 Aq 
4tt 


f A r , dx I p(^) djJ / chrom ( A, n')$(A - A 0 ) dA . 


(I) 
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where h is the Planck constant, B \ 2 the Einstein coefficient for the line, Aq the rest 
value for the central wavelength A of the La transition, and N \ the number density of 
hydrogen atoms in the ground level; the unit vector n is along the line of sight x and 
the unit vector ri is along the direction of the incident radiation; p(^) d J - where 
J is the solid angle around n' - is the probability that a photon travelling along the 
direction n was travelling, before scattering, along the direction n ; ; fi is the solid 
angle subtended by the chromosphere at the point of scattering; / C hrom is the exciting 
chromospheric radiation and T is the coronal absorption profile. In the following 
we assume that the intensity of the chromospheric La radiation is constant across 
the solar disk and that the velocity distribution of the scattering hydrogen atoms 
is Maxwellian. The dependence of the La scattering process on the angle has 
been taken from Beckers and Chipman ( 1974) and we adopted the value given by 
Gabriel (1971) for the ratio between the neutral hydrogen density and the proton 
density at different temperatures (because of the low coronal density all hydrogen 
atoms are assumed to be in the ground level, therefore Ah /X p = X\/X p ). This is 
not entirely correct, since temperatures in our model are 'effective’ temperatures, 
resulting from the polytropic index used in the energy equation. We will come back 
to this point in Section 6. 

In order to evaluate the La intensity in coronal streamers from Equation ( 1 ), we 
need to know how streamers extend in the third dimension. Because our models 
are axisymmetric, it is realistic to assume that they give the distribution of physical 
parameters in a meridional plane, identiried with the plane of the sky and normal to 
the line of sight. If we focus on the dipolar model, we recognize that its geometrical 
configuration is highly reminiscent of the conditions observed at solar minimum, 
when streamers are concentrated along the equator. Therefore, as a first hypothesis, 
we assume that streamers extend all the way around the equator, in a continuous 
belt, and calculate, on the basis of Equation (1), the radial distribution of La 
intensity in a dipolar geometry. To this end, in the following, densities along 
the line of sight are considered equal to those given by the model at the same 
latitude and radial distance. Figure 9 gives the radial profile (solid line) of the La 
intensity, evaluated along the streamer axis (which, in a dipolar geometry, lies in 
the equatorial plane) up to a height of 4.5 /?<.,. Values at larger distances are not 
given, since, beyond that height, field lines are open and the La brightness would 
no longer originate in the streamer. Moreover, open-field regions are not described 
realistically in our simulations, their density being definitely overestimated (see 
Wang et a/., 1992, 1993, for further comments on this point). This is apparent also 
from the slope of the intensity vs distance curve, which keeps constant over all 
the computational domain as if densities decrease linearly with distance. This is 
unexpected, as the line of sight, in regions close to the cusp height and beyond it, 
crosses mostly through the low-density open-field regions. Altogether, the slope of 
the La intensity gradient predicted in the case of a continuous belt of streamers 
circling the Sun, is open to criticism. 

The inaccuracy of the brightness vs distance profile may be ascribed both to 
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Fig. 9 Predicted La intensity vs radial distance along the axis ot an equatorial streamer in a dipolar 
geometry, lor a plasma .i = I 0 (solid linci. The streamer is assumed to extend all the way around 
the equator, and the calculation is performed up to a height of 4 as at higher altitudes most of 
the contribution to La intensities will come from open regions. Results from rocket observations of 
La intensities Irom a quiet coronal region arc also shown (dots). 


an incorrect treatment ot the open field region and to an inappropriate geometry. 
The latter factor would be influential if only the slope at large distances is wrong. 
However, the values of the La intensity at low heights also look too high, in 
comparison with the tew data points available so far. In Figure 9 we show (dots) the 
La intensities given in Figure 6-1 of the UVCS Science Requirement Document: 
these data refer to a quiet region and were obtained in 1979 during a rocket flight 
(Kohl et nl., 1980). Because they do not refer to a streamer, they do not provide 
any information about the change in the slope of the La intensity gradient in the 
cusp region, but we can infer that La streamers would be a few times brighter than 
appears from the rocket data. The La brightness predicted from the model is about 
one order of magnitude larger than observed in quiet regions, which, taking into 
account that our densities are realistic, is too large a factor. This rules out our initial 
hypothesis of a continuous belt of streamers circling the solar equator. It is worth 
pointing out that our conclusion is consistent with observations, as these seem to 
indicate the presence of several streamers spread, at a given latitude, over different 
longitudes (Dollfus, private communication). 

From this analysis we conclude that we need both a better simulation of open 
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tield regions and a definition of the streamer geometry through a full 3-D model 
However, we may still get realistic predictions by adopting a priori the streamer 
configuration and neglecting the contribution from the outer low-density regions. 
The resulting profiles of intensity vs distance will illustrate the behavior of the 
La brightness for different geometries and UVCS data will eventually allow us to 
identity the more realistic configurations. 

Figure 10 shows the La intensity vs distance profiles assuming three different 
streamer configurations in a direction normal to the streamer axis and parallel to 
the line of sight (we remind the reader that streamers are supposed to lie in a 
meridional plane). Either the unqulur width of a streamer is constant with height 
(fan-shaped streamer) and equal to its base angular width (as seen in the plane of 
the sky) or its ! in tar width is constant with height (constant-thickness streamer) 
and equal to the distance between the footpoints of the highest closed field lines 
(as seen in the plane of the sky) or streamers are cone-shaped’ structures and 
therefore have a width, initially equal to that of a constant-thickness streamer, 
which decreases linearly with height up to the cone vertex, identified with the 
streamer s cusp. Figure It) shows that, in a dipolar geometry, an equatorial fan- 
shaped streamer differs negligibly from a continuous belt of streamers circling the 
Sun. This is due to the large width of the streamer (half width % 38°): outer regions 
contribute to the emergent intensities only at great distances where densities are 
too low to affect significantly the La brightness. As a consequence, unless future 
observations will show streamers to be brighter than assumed so far, we are led to 
discard the hypothesis of streamers as constant angular width structures. On the 
contrary, constant thickness and cone-shaped streamer structures lie close enough 
to the observed data points to be equally plausible. 

However, our model locates the cusp only approximately, both because our 
model does not take diffusive effects into account and because of the coarse 
resolution of our mesh points. Hence, in order to illustrate, in a cone-shaped 
geometry, how different cusp altitudes affect the La intensity gradient, we have 
considered the cusp height as a free parameter and evaluated the resulting radial 
profiles in the usual dipolar geometry and high /3 plasma. Figure 1 1 shows the La 
intensity vs distance profiles for a cone-shaped streamer whose vertex - i.e., cusp 
height - is located at altitudes ranging between 2.5 and 6 R !q. We point out that such 
a large variation in the position of the streamer’s cusp far exceeds the uncertainty 
of the model and is shown only for display purposes (although the procedure is not 
entirely consistent, as different J values would be required to build models with 
such different cusp heights). We conclude that the La intensity and the slope of 
the La intensity gradient initially (i.e.. close to the Sun) depend only weakly on 
the shape of the streamer, but, at large distances, are dictated by the streamer’s 3-D 
structure and. in the case of a cone-shaped feature, by its cusp height. 

So far, our examples referred to a global dipolar streamer. However, our results 
can be extended to the quadrupolar and hexapolar model configurations by taking 
into account the differences in the streamer geometry. As we have shown in Sec- 
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Fig. 1 0. Predicted La intensity vs radial distance aiong the axis of an equatorial streamer in a dipolar 
geometry, tor a plasma J = 1.0. The streamer thickness in a direction normal to the plane of the 
sky is assumed either to have a constant angular width (fan-shaped streamer) equal to its aneular 
base width (^76 ) or to be constant with height (constant thickness streamer) and equal to its base 
width (as I /?:,) or to decrease linearly with height (cone-shaped streamer) up to the cusp height 


tion 3, densities are about equal tn low and high latitude streamers. Therefore, the 
La brightness from these features turns out to scale, with respect to that originat- 
ing in a dipolar equatorial configuration, in the same proportion as the streamer 
thickness. The La brightness from an equatorial hexapolar streamer, for instance, 
will be a factor 2—3 smaller than that from an equatorial dipolar streamer and about 
equal to the brightness from the high-latitude hexapolar streamer. 

We did not consider, yet, the case of off-axis observations, which should provide 
a more comprehensive test of the model by allowing us to determine the physical 
parameters of the streamer across its axis, ov er a meridional plane. For instance, if 
SOHO UVCS were to observe an equatorial streamer - symmetrical with respect 
to a meridional plane through its axis - when its symmetry plane lies in the plane of 
the sky, it should be possible, via off-equator observations, to check the shape and 
physical parameters of the streamer in the meridional plane purportedly described 
by the model. 

Figure 1 2 gives the La intensity gradient, measured in the plane of the sky along 
directions parallel to the axis of the streamer, lor the usual dipolar configuration 
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(0 1 .0, /? cusp _ 4.5 R : ), in the hypothesis of a cone-shaped structure. The La 

intensity vs distance gradient is here dictated by the variation of density across the 
streamer and by the decrease of the integration length along the line of sight as 
we move off the equatorial plane. The latter factor is responsible for the intensity 
drop which occurs at increasingly lower distances, as we move to larger offsets, 

and mimics the behavior observed in the cusp region in the equatorial plane also 
shown in the figure. 

If we now move in the plane of the sky, along a direction normal to the streamer 
axis, and we evaluate the La intensity at increasing offsets, we would guess 
t at the La brightness decreases proportionally to the decrease in the integration 
length, and that, as a consequence, the ratio of La intensities evaluated at positions 
corresponding to increasing offsets can never be lower than the ratio between the 
corresponding integration lengths. For instance, at a distance - measured along the 
axjs of the streamer - of 1.25 f?, ? , the intensity ratio I R = I Lo 0 2 % // Lq 0 4 R 
is s=l.5, while the ratio between the integration lengths at those offsets is 2=1.3. 
The decrease of density, as we move off axis, accounts for a 10% increase in the 
intensity ratio over the value predicted on the basis of the ratio between integration 
lengths. However, the intensity ratio is smaller than the ratio between integration 
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Fig. 12. Predicted La intensity vs distance (measured along the axis of an equatorial streamer) for 
different oil-equator ottsets in a cone-shaped equatorial streamer: solid line, intensity along the axis 
of the streamer, zero offset: A. intensity for an offset of 0.2 fl: : x, intensity for an offset of 0.4 R 
As usual the conliguration is dipolar, the base width is %38° and the cusp height is 4.5 R 



lengths whenever the shorter integration path, at a larger offset, runs either through 
the density enhancement at the edge of the streamer (see Figure 5), or through 
regions with a lower temperature (and, as a consequence, a higher percentage of 
neutral hydrogen atoms). At a distance of ^2 R . , this effect makes the intensity 
ratio about 10% smaller than the ratio between integration lengths. Thus, from 
Figure 1 2, we conclude that the ratio between La intensities, evaluated at the same 
distance (along the streamer’s axis) and different offsets, is approximately equal 
to the ratio of the integration lengths. In the hypothesis of an axially-symmetnc 
density distribution, we conclude that densities play a secondary role, with respect 
to the streamer s geometry, in determining the value ot La off -equator intensities. 
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5. La Emission from a Rotating Streamer 


The extended lifetime of SOHO makes it possible to observe streamers over a 
prolonged period of time as they are carried around by solar rotation Hence, 
generally, the line of sight will cut through the structure obliquely and the angle 
between the streamer's axis and the integration path will vary from 0° to 90 as the 
streamer moves from the central to the limb meridian. When the line of sight cuts 
normally through the streamer, the region which is closest to the Sun and, therefore, 
has the highest density, is located at the streamer axis: as a consequence, the highest 
contribution to the emergent intensity comes from this region. On the contrary, when 
the line of sight cuts obliquely through the structure it may happen, depending on 
the angle between the line of sight and the streamer axis and on the axial vs 
transverse density gradient, that the highest contribution to the emergent intensity 
originates from a region at some distance from the streamer axis, //evolutionary 
effects are negligible - i.e„ if streamers are stable throughout a period of time 

- we may use this effect to get information on the density prohle in a direction 
normal to their axis and. through prolonged observations, eventually reconstitute 
their entire structure. In other words, for stable structures solar rotation allows us 
to see streamers under different perspectives and use tomographic techniques to 
obtain their 3-D configurations. 

In the previous section, we assumed that the streamers footpoints were rooted 
at the same longitude. However, if streamers are rooted in active regions, it is 
likely that their footpoints are rooted at the same latitude - say, along the equator 

- inasmuch as positive and negative polant.es tend to align along the east-west 
direction Although our model seems inappropriate to deal with this case since 
it is not realistic to have magnetic ‘poles’ along the solar equator - as long as 
we do not have a 3-D simulation it is plausible to focus on the streamer sector 
and adopt the representation provided by the model to describe streamers tying on 
the equatorial plane. This allows us to explore the capabilities of the tomographic 
technique, because, in this hypothesis, the model provides a complete description 
of the behavior of density along the line of sight (at least for on-axis observations). 
Hence, in the following, contrary to what has been hypothesized so far, the streamer 

is assumed to lie on the equatorial plane. to 

Figure 13 shows how individual elements along the line of sigh c 
the total La intensity measured in the equatorial plane at a distance of 1.25 K®. 
When the streamer is at the limb (streamer longitude 90°), its axis lies in the phme 
of the sky and is perpendicular to the line of sight. Therefore, the element lying 
at l 95 flc, along the axis is the element closest to the observer and provides the 
highest contribution to the emergent intensity (top left panel of Figure 13). As the 
streamer is carried around by solar rotation, different elements, at some distance 
form the axis, become the major contributors to the total intensity^ ig 
demonstrates the progressive shift of the element which most 
emergent La intensity, as the streamer longitude changes by 30 . The bottom right 
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panel, for instance, shows that, when the streamer has rotated by 30° behind the 
plane of the sky, the highest contribution to the emergent intensity comes from the 
element at a distance of 0.72 R , along the line of sight (distances along the line 
of sight are counted from the streamer axis). Thus, by taking an extended set of 
data, at different locations along the axis of the streamer and at different rotation 
angles, we eventually get a complete map of the density of the structure. 

We do not give any further example of this technique as a more realistic choice 
of boundary conditions in the open-field regions will modify the distributions of 
Figure 13. We anticipate that a lower density in open-field areas will result in a 
steeper decline of the contribution from elements located outside, or at the outer 
edge of the streamer. This effect may help getting a density map with higher 
spatial resolution than otherwise possible. Although stable structures may be a 
minority within the streamer family, the example of Figure 1 3 shows that it is worth 
developing this methodology further, as a means for an observational determination 
of the 3-D streamer's structure. 



6. Discussion and Conclusions 

Our purpose has been twofold: to provide guidelines for UVCS observational 
sequences and to compare our simulations with the scanty data available in order 
to guide further development of the numerical model toward more realistic global 
configurations. 

The first objective has been reached, insofar as we have presented a set of 
predicted profiles of La intensity vs distance, both for on-axis and for off-axis 
observations and for different rotation angles of the streamer. These profiles can 
be easily adapted to different magnetic configurations. On the basis of our model 
for an axially-symmetric structure, we also show that off-axis observations allow 
an identification of the streamer’s dimension along the line of sight. Finally, we 
have shown that this capability, combined with prolonged observations of a stable 
streamer at different longitudes, leads to a 3-D map of densities in streamers for 
comparison with our global simulation. 

The model uses a polytropic relationship between density and pressure, rather 
than a full energy equation. Hence, the temperatures we predict are effective tem- 
peratures. Observationally, the La brightness depends on the electron tempera- 
ture, via the neutral hydrogen abundance, and on the kinetic temperature, via the 
coronal absorption profile. As long as we consider integrated La brightness, the 
effect of an incorrect absorption profile is probably negligible. Model tempera- 
tures are, however, lower than temperatures derived from streamer observations 
(Liebenberg, Bessey, and Watson, 1975), so we apparently overestimate the neutral 
hydrogen abundance. Nevertheless, measured temperatures in streamers have such 
great uncertainty that we cannot resolve this issue until UVCS provides accurate 
measurements of the electron temperature. 

Ultimately, our second objective is the more relevant. The simulations point 
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to the need for (a) a different choice of boundary conditions in open regions than 
inside the streamer, (b) an extension of the model to different 6 values, and (c) the 
development of a 3-D model. We are presently working on these issues and expect 
to get a realistic simulation of open held regions shortly. Our goal is to attain 
a reasonable global model that simulates both streamers and coronal holes by 
achieving agreement between model predictions and existing observations. This is 
a method to fully determine the 3-D structure of streamers to within the resolution 
of the model and observations. 

Traditionally, three-dimensional mapping of the electron density in coronal 
structures has been performed via the image reconstruction technique originally 
developed for X-ray tomography (Altschuler and Perry. 1 972; Perry and Altschuler, 
1973; Wilson, 1977), or via the already mentioned Van de Hulst’s ( 1 950) method. 
Both procedures have been subject to criticisms (Bagenal and Gibson, 1991) be- 
cause of their heavy computational requirements. While the ideal method for a 3-D 
mapping has yet to be found, an alternative approach to devise theoretical models 
with free parameters that are calibrated against observational data, is being devel- 
oped (Bagenal and Gibson. 1 99 1 ). The present work uses this alternative approach 
and the results achieved so far indicate that it is worth pursuing. We conclude 
by pointing out that an agreement between model predictions and observations 
will allow us to determine the magnetic field vector throughout the streamer, thus 
providing a complete picture of these so far elusive structures. 
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Abstract. Wc describe a two-dimensional time-dependent. numerical, magnetohydrodynamic model 
for the determination o I the physical properties of coronal streamers from the top of the transition 
/one [R = l) to 15 R. . Four examples are given; for dipole, quadrupole. and hexapole initial 
Held topologies I he computed parameters are density, temperature, velocity, and magnetic field. In 
addition to the properties of the solutions, their accuracy is discussed. We use the model as the basis 
ior a general discussion of the way boundary conditions are specified in this and similar simulations. 

1. Introduction 

We present results from a recently-developed numerical model of coronal struc- 
ture. The immediate reasons for a new model were to extend the outer boundary 
farther from the Sun and to gain the experience necessary for development of a 
three-dimensional model. A result of this process has been a close examination of 
the physical details of the solution and how they depend on the way the boundary 
conditions are specified, An immediate application will be the simulation of stream- 
ers in support of the Ultraviolet Coronagraph and Spectroheliograph (UVCS) and 
the Large Angle Spectrometric Coronagraph (LASCO) on the Solar Heliospheric 
Observatory (SOHO). These instruments will be able to measure the temperature, 
density, and flow vector in the corona. With model calculations, it will be possible, 
for example, to estimate the magnetic field vector. 

Numerical models of coronal structure have been published sporadically, at 
long intervals, over the past twenty years. The first (Pneuman and Kopp, 1971) 
demonstrated the feasibility of such models, treating isothermal flow and arriving 
at the solution by iterating on the electrical currents. However, a more efficient 
and flexible method is to consider an initial-boundary value problem in which 
the steady state is found holding the boundary conditions constant and allowing 
the solution to relax in time from an essentially arbitrary initial state. Steinolfson, 
Suess, and Wu (1982) applied this later technique to the analysis of a polytropic 
dipole configuration for a range of plasma ,j (ratio of internal pressure to mag- 
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netic pressure). Sleinolfson (1989. 1991 ) and Guo et al. (1992) have used this 
steady-state solution as the basts tor studying coronal mass ejections and streamer 
evolution with shear, which can be simulated using a nearly identical numerical 
model. Details of the numerical schemes and results can be found in the referenced 
publications. 

We revisited this problem for the reasons mentioned above. However, we also 
consider that such complex numerical models are rarely without problems or uncer- 
tainties. When the models are used for analysis of data and for predictions, the only 
reliable validation is to develop an independent model and compare the results. 
E\en when both for all) models are fundamentally correct, this process generally 
leads to new or deeper understanding of the problem. In the present case, this is 
precisely what has happened. We have gamed a better insight into the physical 
basis of the criteria which should be adopted in specifying boundary conditions. 
The results from this constitute an important part of the present study. 

The physical and numerical simulation is described in Section 2. Section 3 
details numerical models of dipole, quadrupole, and hexapole magnetic fields. 
Section 4 is a discussion of numerical precision of the solution and the boundary 
conditions, putting the discussion into context with earlier models so far as is 
possible. Section 5 contains our summary and conclusions. 




2. The Physical and Numerical Simulation 

We assume axisymmetric, single fluid, polytropic, time-dependent ideal magneto- 
hydrodynamic flow and perform the calculation in a meridional plane defined by 
the rotational symmetry axis of the magnetic field. The coordinates are (r, 9 , <p) 
with 0 being the ignorable coordinate. For the magnetic field boundary condition, 
we take the radial lield component at the lower boundary to be that given by a 
vacuum dipole, quadrupole, or hexapole potential magnetic field. The flow there- 
fore has reflective symmetry across the equator and the calculation need be done 
in only one quadrant. The equations of motion that describe this flow are: 
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The dependent variables are the density, p. the pressure, p. the radial and 
meridional velocities. r r and r„, and the radial and meridional magnetic helds 
B r and Bo. The constants M. . CV. ,. and p are the solar mass, gravitational 
constant, the poJvtropic index and the magnetic permeability. 

These equations are solved in a computational domain extending from the Sun 
( I R : ) to 1 5 R . . from the pole to the equator. It is assumed that meridional flow 
is zero at the pole and equator. The grid is divided so that there are 37 endpoints 
in the radial direction and 22 endpoints in the meridional direction, witluhe radial 
grid size slowly increasing with radius. The meridional grid is divided so that 
points lie equidistant on either side of 0 = 0 and 0 = 90°, at p _ 050 

2.25°, 6.75°, . ., 87.75°. 92.25° The algorithm adopted here is the Full-Implicit 
Continuous Eulerian (FICE) scheme described by Hu and Wu (1984). For time 
stepping a second-order accurate forward differencing scheme is used, with the 
step size being of the same order as given by the Courant condition because the 
magnetic field is calculated explicitly. Smoothing is used when gradients become 
too large, i.e.. at shocks (which do not occur here). At the inner boundary, the 
flow is subsonic and sub-Alfvenic so that two of the six independent variables are 
calculated using compatibility relations iHu and Wu, 1984). A brief summary of 
the compatibility conditions for the present model is given in the Appendix, along 
with details on how the boundary values and conditions are applied. We choose 
to specify the radial and meridional magnetic fields, temperature, and density. The 
radial and meridional How speeds are computed from compatibility relations (i.e.. 
Equations (A. 1 ) and (A. 2)). At the outer boundary, the flow is restricted to being 
both supersonic and super-Alfvenic. In this case, all variables at that boundary can 
be calculated by simple linear extrapolation from the first (or first two) grid points 
inside the boundary. In this study, we did not perform the comparison between 
the present boundary conditions and conventional boundary conditions. However, 
in a recent study by Sun (1991 ). it was shown that the statement of the boundary 
conditions in the Appendix eliminates the spurious waves generated by boundary 
disturbances and which can cause numerical instability. 

We start with an essentially arbitrary initial state and allow the flow to relax 
in time while holding the boundary values constant. In the present case the initial 
flow field is a polytropic, hydrodynamic solution to the steady-state radial flow 
equation of motion (e.g., Parker, 1963) Miperimposed on a potential magnetic field. 
That this is neither a self-consistent nor stable solution to the steady-state MHD 
equations is irrelevant since the How is allowed to evolve in time under the control 
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Ot the equations ot motion. The main concerns are that the numerical solution 
be stable and ot sufficient accuracy to define the physically interesting aspects of 
the solution, and that the relaxation proceed long enough that an acceptably close 

approximation to the steady state has been reached. We address these issues briefly 
in Section 4. 



3. Detailed Results from Four Specific Models 


We report here on four specific models. The results are grouped hrst accordins 
to the way in which the physical variables are plotted li.e.. cither versus radius 
or versus polar angle) and second according to which of the four examples the 
plot is for. In these tour examples, three magnetic held geometries are used: a 
dipole, a quadru pole, and a hexapole; the scalar potentials are therefore proportional 
to P> leostf). P\ (cosd). and Pj, Icoshi. respectively, where P„ (cos#) is the 
Legendre polynomial of degree n. There are two dimensionless free parameters: 
the polytropic index. - , and the plasma i. We use = 1.05 in all cases. J = 1,0 
lor all three held geometries, and. in addition, do a dipole calculation for A = 0 
In these case, i is evaluated at 1.0 /?. at the equator, where the held strenmhTs 
1.67 G both for ,1 = 1 .0 and .1 = 0.2. For the high i cases, the base temperature 
and density are 1.8 x 10'- K and 2.25 • It)-'' cm \ For the low .j case, they are 
1.44 x 10 ’ K and 5.61 10 ' cm The three magnetic held geometries naturally 

lead to a single equatorial streamer, a mid-latitude streamer, and both an equatorial 
and a mid-latitude streamer for the dipole, quadrupole and hexapole. respectively. 

Results from the four examples w ill be referred to as follows: 

(a) Dipole, A — I (). 

(b) Quadrupole, A = 1.0. 

(c) Flexapole, A = 1 .0. 

(d) Dipole. A = 0.2, 


The initial state temperature, density, and velocity profiles are shown in Figure 1 
The temperature curves appear irregular due to the small change in temperature 
over the relatively large radial range - a consequence of the polvtropic index being 
near unity. Only three significant figures were retained after the calculation so what 

is seen here is roundoff error in the plotted results rather than in the computed 
results. 

The hnal, steady-state magnetic held geometries for the four cases are shown in 
Figure 2. Here is seen the well-known property that the flow is nearly radial beyond 
3-4 R,~„ The flow is field-aligned everywhere and field lines which cross the outer 
boundary reach to oo. The streamers are those volumes which are magnetically 
closed fi.e., the field lines return to the surface of the Sun) and it is evident that 
relatively small volumes in the streamers remain magnetically closed in comparison 
to the initial state where all Held lines were closed. These closed volumes are 
surrounded by a low density shell but. as will be shown below, the densities in the 
large coronal hole-like open regions are otherwise only slightly lower than in the 
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L-0 Initial State, M,o ‘“0 Initial Slate, ^0Z 




Fig. I . Density, temperature, and velocity profiles in radius that were used for the initial (t = ()) 
state in the relaxation. On the left arc the protiles tor the .i = 1.0 cases and on the rieht arc the 
profiles tor the J = 0.2 case. Note that, except lor the velocity scales, the scales differ between the 
two panels. Because the polytropic index is near unity, temperature changes slowly with radius and 
the irregularities in the temperature profiles should he interpreted as noise. 


streamers. In each panel of Figure 2. four dashed lines are shown and labelled A. 
B, C, or D. These lines indicate the radial directions used below to plot variables 
versus radius. 

The physical times allowed for the relaxation in these four examples were: 
(a) 22.22 hours for the :j = 1 .0 dipole; (b) 1 6.67 hours for the li = 1.0 quadruple; 
(c) 1 8.06 hours for the ii = 1 .0 hexapole; (d) 1 9.44 hours for the i3 = 0.2 dipole. 
These times are determined by how long it takes for any fluctuation to be advected 
out through the outer boundary of the solution domain. This in turn depends on how 
large the flow speed is and whether the fluctuations represent inward propagating 
waves. In general, the times listed above are the minimum required for a stationary 
fluctuation (i.e., non-propagating in the solar wind frame) to be advected from 
I R , to 1 5 R~. at a typical flow speed in the open regions. This sometimes leads to 
small residuals in the relaxation near the outer boundary at 15 fto, but the solutions 
inside 7 /? r, that are shown here are quite steady. This is another point that will be 
reviewed in Section 4. 

Figures 3 and 4 are plots of density and radial velocity versus radius. The plots 
are made in the directions indicated in Figure 2 so that, for example, in each panel 
of Figure 3 the density is plotted in the four directions A, B, C, and D indicated 
in the corresponding panel of Figure 2. In both of Figures 3 and 4, the four panels 
corresponding to the four panels in Figure 2 are clearly labeled. The density profiles 
have been divided by their corresponding initial state (t = 0) profiles from Figure 1 
because the density changes by several orders of magnitude between the Sun and 
15 Ri7y . The plots here extend only to 7 /?, r . because there is no new information 
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h'hplOH Th CnS '^ aS 3 * unctlon °* rac ^ us - panel is lor the corresponding case in Figure 2 as 
! bnhe J L | C n h CS ^ ^.Tu a ! 0ng lhe dirca,ons sh °w„ F^re 2 For example, the tour cuAes 
ot Fieure 7 n h 6 J ^ c cd A ' B ' C - D ' arc alon g «hc lour directions shown in the third panel 

Film re i 1 “* ? labelled ln Ihe same manner. Each curve has been divided by the initial profile (see 
ien«,v ‘ y enhanccmcnl ls ' nd ' ca '^ by values greater than unity, and mce versa The 

50% abo^eTSal state S ‘ reamerS ^ C ' Car ' y V,S,We - gCnera " y bein ® on lhe orde ^ of 25% to 


contained outside this radius - the flow is already supersonic and essentially radial. 

umtng briefly to each figure individually, we begin by noting that a density 
enhancement is indicated by values greater than unitv. and vice versa. The density 

oHWtnsZ 5 S ‘ rCamer ! hCre 3re Clearly vis,ble - generally being on the order 
ot 25% to 50% above the in.Ual state. The base density for the d = 1 .0 cases is close 

to that reported by Allen (1955) for the base of the quiet corona and the dens.ty 
prohle shown here has generally the r.ght behavior for streamers - as shown by 

TT- T™* (a) ’ (C) ’ and (d) ’ and curve D in case (b) - Curve D for ca se (a), 
e • .0 dipole, is an example of the density deficit on the flank of a streamer 

a is typical of the results for all the examples. In contrast, the density in the 
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A 





12345G70 1 2 3 4 5 0 7 0 

Solar Radii Solar Radii 

Fig. 4. Radial velocity as a Junction ot radius. Each panel is tor the corresponding case in Figure 2. 
The curves are plotted along the directions shown in Figure 2, as in Figure 3. The velocity inside the 
streamers is seen to be essentially zero. 


centers of the open regions (curve B in all cases, curve C in case (b), and curve D 
in case (c)) is little different from the initial state, being only slightly smaller. This 
is only surprising when comparison is made to coronal hole observations (Munro 
and Jackson, 1977) wherein the density was reported to be more than an order of 
magnitude less than m streamers. This difference is a natural consequence of the 
properties of a polytropic model and the choice we have made for the boundary 
conditions on temperature and density - that they be independent of polar angle. 
The choice leads to both the high density shown here and the low flow speeds 
shown below on open field lines, irrespective of the open streamline geometry. 
To model true coronal hole flow with a polytropic gas would require at least an 
elevated temperature in the open regions and probably also a lower density at the 
base (Suess et a/., 1 977; Suess, 1 979). 

The radial velocity is shown in Figure 4, at the positions indicated in Figure 2. 
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s A peed rihf, ab0Ve ' Md “ ,S ge " erally "" “ se ,n Polytrop,c models. ,he Now 
m Figure L Ke Tamer T f “ ,hs ; ndis " ,rt >' d «»» speed shown 

above „ e d ° *» lh = «P« *« Teg, on 

this - rnH rhf ft “ ?• he J - 0 2 d ‘P° le ,s ‘he most extreme example of 

v^es 0 bv 0 n„ P |v , , rt ! e KmPera,Ur t S '" Ce - due to ,he Poly*"** being LOS. 
this is an >iTa V ^ PerCent lhrou 3 bour the compulation domain. However 

dissipation mechanisms ^^Howlver.^has^bc^sbow^th^^ a^pol^ropic^ndex 
on the order 1 .1,5 is reared reproduce observations of comnalden^eMSue" 

How!"dh™' ma r e,,Cally ° pcn reg,ons - J bhough euqivalen, to coronal hole 
> simulate coronal holes because the How speeds are far too small 

the ? b " reaso " ilble liow speeds in this model it would be necessary to have 
he temperature vary across the base of the open region - which is well within 
he capability of the model. Such a variation has bL shown ,o leZduce all 

geom™rTw,Ihmeeff 0 ', COr ° nal h °' e "° W ^ '“ d 10 accuraK simulations of the 

at the ed^eTsuess^ ‘ '*"* ° f ,he hole tha " 

I . & ^ ' n con trast to the open regions, the densities in the 

regions are similar to observed streamerdensities and we feel this model is 
dterefore a good upproximalion to streamer geometry. The lempem urn mu « s , 
be qualified as an effective temperature, but can be used for diag^c puZses n 
combination with planned observations on SOHO/UVCS ^ 

nlnrtT ° f '** ^ be be,ter V,ewed and ™re easily understood when 
plotted versus polar angle at different heliocentric distances, than versus radius at 

b H P a " EleS Such pl0B are show " '°T the density, radial veToaty and 
total held strength in Figures 5, 6, and 7, respect, velv. * 

the ^ - C iVh° W i th t de " SUy dr ° P ad J acent to the streamer. In the panel for 
dTni ~ ° d,p ° e ’ thlS dr ° P is quite lar S e - we| l resolved, and iead S P into the 

occur fsa" the^e"' ^ ^ Uat ° rial streamer The only place this does not 
p _. h 6 u “ Where the dens,tv ,s held constant. The width of the density 

enhemcement ,n the streamer decreases with height, just as the width of the streamer 

hself decreases with height (e.g„ Figure 2). Essent.ally the same thing is seen for the 

higher Ldwrrtn^ < q H Uant,tatl ve differences: (,) The streamer is much 

.gher and w.der. („) The density depletion on the Hanks has a smaller amplitude 
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Fig. 5. Density versus polar angle, between the pole (0° ) and the equator (90°). Each of the curves 
is labelled according to the heliocentric distance it refers to. Thus, the curves labelled 1.70 ft indicate 
the density at 1 .70 ft,: heliocentric radius. The density at the base is constant and so the curves there 
are flat. Above the base, there is a small density enhancement in the streamer (ca. 5% to 50%) and 
a trough in density at the edge ot the streamer. In the middle of the open region, the density is very 
close to what it was in the initial state (see also Figure 3) The reason it is not small is that we have 
used constant temperature and density at the base. To produce a true coronal hole-like profile would 
have required at least an increase in the temperature at the base of the open region (Suess et al. 
1977). 


These differences are the primary reason we conclude that solar streamers are better 
described by a 3 ~ 1.0 plasma than by a A = 0.2 plasma. Qualitatively, a similar 
result is found for the quadrupoie and hexapole. However, it is obvious that the 
hexapole is only marginally resolved with the present grid density - there is really 
only one meridional grid point inside the mid-latitude streamer at any given height. 
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The radial veloc.ty in Figure 6 drops prec.p.touslv from the magnetically open 
reg.on to the inside of the streamer. That the velocity is not identically zero' ms.de 
the streamer is a result of numerical diffus.on and is a measure of this numerical 
artifact in the PICE scheme. For example, at 1 .70 ft, , n the .3 = 1 .0 dipl the 
\e ocity drops from ca. 60 km s 1 at the edge of the streamer to about 3 km s“‘ 
.ns.de the streamer. 3 km s^ 1 is hardly above the noise level in the plots and the 
associated kinetic energy is too small to affect the dynamics of the solution Such 
slippage will, nevertheless, occur in all numerical solutions. At larger heights 
le.g 4.90 and 7.14 ft,) there is small, but finite flow near and in the neutral sheet 
dividing regions of opposite magnetic polarity. This is qualitatively like what is 
observed in the solar wind in the interplanetary medium. The .3 = 0.2 dipole again 
exhibits properties unlike the Sun in the sense that the verv low tiow speeds inside 

the streamer seem to still exist even at 7. 14 ft . - far outside the observed extent 
of closed streamers. 

Figure 7 shows the variation of the total magnetic field strength, l B; + B -) 1 -- 
across the streamers. The most interesting thing to note in these plots is the enhance- 
ment in total field strength on the Hanks of the streamers. This is what -confines’ 
the streamers. The held strength for the .3 = 0.2 dipole is seen to vary smoothly 
with little distinct evidence of the streamer. This is just another indication that the 
presence of the plasma has had little effect on the held geometry in this low-3 case. 

4. Accuracy and Stability of Calculations 

This numerical model has been found to be weakly subject to the Courant condition 
on size of time step. Therefore, the size of the time step decreases as the largest 
values of the temperature and magnetic held increase - along with the maximum 
sound and Alfven speeds anywhere in the grid. Counteracting this, the higher 
characteristic speeds lead to a somewhat faster relaxation time. However, generally 
shorter time steps are required for smaller 3 calculations. The flow speed also 
plays an important role in determining the relaxation time to a steady state - the 
initial state is a disequilibrium conhguration. This imbalance must have time to be 
advected from the base through the outer boundary. The physical time this takes can 
be estimated by taking a typical (but small) value for the flow speed and caculating 
how long it would take the plasma to flow at this speed from the base to the outer 
boundary. For example, at 150 km s' 1 , to 15 ft. , this takes 18 hours (relaxation 
times we have used here are given in Figure 1 ). 

A second consideration is gndpoint resolution. The grid used in these examples 
is 4.5 in latitude and about 0.24 ft 0 in radius near the base - increasing slowly 
with radius. This is sufficient to adequately resolve the geometry and flow on the 
scale shown in Figure 2. However, if finer scale information is required in, for 
example, the core of the streamers, a denser grid would be required. 

Always a serious consideration in these time-dependent. non-Cartesian MHD 
calculations is the conservation of magnetic flux - that V • B = 0 is maintained 
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Fig. 6. Radial velocity versus polar angle, between the pole and the equator. Bach curve, plotlcd 
lor different heliocentric distance, is labelled in the same manner as in Figure 5. The velocity in 
the magnetically closed regions is essentially zero. The reason it is not identically zero is that there 
is a smaU amount ot numerical diffusion - quite small as indicated by the velocity being less than 
10 km s inside the ti = 10 dipole streamer at 2.30 



Fig, 7. The to: 
Each curve, p U 
The held in th 
amplitude, as \ 



at all times. The condition is maintained here through accurate differencing rather 
than a self-correcting scheme. No anomalous acceleration due to errors in flux 
conservation is apparent in the results. The numerical scheme is pressure-based so 
it is limited by stability to large and moderate .j values (e.g., J > 0. 1) - which 
turns out to be the same restriction for maintaining V ■ B = 0 to the required 
degree. 
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DiDole.3eta= 1.0 
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Hexapole.3eta = ] .0 
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Fig. 7. The total magnetic field. ( /i; + Bi)'*-. versus polar angle, between the pole and the equator 
Each curve, plotted tor different heliocentric distance, is labelled in the same manner as in Figure 5 
The field in the v.cinity of the current sheet above the cusp in the streamers has a greatly reduced 
amplitude, as would be expected. The effect is amplified above the mid-latitude streamers.' 
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reduces to v • V(p/p 7 ) = 0 when a steady state is reached, which means that ( p /pi ) 
is then a streamline constant. This becomes an analytic test of the achievement of a 
steady-state solution in our case. The boundary values of p and p are the same at all 
latitudes. Therefore, (p/ p 7 ) = 0 has the same value everywhere in the computation 
regime as it has on the boundary if a steady state has been reached. We have checked 
this for the cases shown in Figure 2 and find that for the dipole and quadrupole it is 
constant to within a maximum of \% and for the hexapole it is constant to within a 
maximum of 4% (average values over the whole grid are less than 1 % in all cases). 
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5. Discussion 

The new feature of this model, with respect to analogous simulations, is the exten- 
sion of the outer boundary to 15/?. This is not a conceptual advance, but this and 
the stability and ruggedness ot the code make it very useful for simulating realistic 
coronal conditions. We present new' results for quadrupole and hexapole fields, 
with their accompanying mid-latitude streamers and open magnetic held regions. 
The Alfven speed ranged between 800 km s and a few tens of km s' 1 . This is 
lower than is believed appropriate for the corona ( Suess. 1 988), but we expect our 
model will now enable simulations with higher Alfven speeds. 

When comparing our results to those of Stemolfson, Suess, and Wu (1982: 
henceforth referred to as SSW). an interesting and important difference becomes 
apparent. In the present calculation, we have held the density and temperature 
constant at the base, allowing the velocity (and, hence, the mass flux) to ‘float’ 
with time in accordance with the compatbility relations determining the velocity 
from the solution inside the computational domain. In contrast, SSW hold the 
temperature and velocity constant at the base and allow the density to change 
according to the compatibility relationships. SSW determine the location of the 
streamer by locating closed field lines and allow ing the velocity to decrease to zero 
at the feet ot these field lines. A consequence is that inside the streamer, the final 
density is considerably higher than the initial density and this is the primary reason 
lor the quantitative differences between their results and ours. 

There is an important consequence of this difference in boundary conditions 
between SSW and the present calculation: the plasma 3 is computed using the 
temperature, density, and magnetic held at the equator and at 1 R^. This is invariant 
in the present calculation, but in SSW this number is different in the final, steady 
state than at the beginning: there 3 was computed using the initial values. Therefore, 
in SSW in the steady-state solution is actually larger than stated for each example 
they did. Thus, our calculation for a dipole with 3 - 0.2 (case (d)) corresponds to 
cases for 3 < 0.1 in SSW. We feel that the way we have done the analysis more 
closely corresponds to what occurs and what is physically known for the Sun and 
therefore leads to a more precise definition of the problem. So, we conclude that the 
present study has demonstrated a preferable treatment of the boundary conditions 
in comparison to earlier calculations. 

A consequence of the precise examples we have done in cases (a) through 
(d), with constant temperature and density, is the flow speed and high density in 
the magnetically open regions - in comparison to what is believed to be the case 
in solar coronal holes. This is a natural consequence of using a polytropic gas 
in which the flow speed is strongly dependent on base temperature. It also does 
not reflect suggestions from analysis of Skylab data that densities at the base of 
coronal holes may be a factor of two smaller than at the base of streamers (G. Noci, 
private communication). In a continuation of this study, we will produce models 
with varying temperature and density at the base. The variation in temperature 
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Will, because it is an 'effective temperature . retiect a difference in energy balance 
and distribution between the base of coronal holes and streamers instead of a true 
temperature difference. 
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Appendix. The Inner Boundary Conditions According to the Projected 
Normal Characteristic Method: a 2D Case 

The inner boundary conditions are obtained according to the method of projected 
characteristics (Nakagawa. Hu, and Wu, 1987) with the FICE algorithm (Hu and 
Wu, 1 984). For the two-dimensional case, the Alfvenic mode does not exist, there- 
fore, there are six eigenvalues. These six eigenvalues lead to six projected normal 
characteristics and to six compatibiliitv equations (see Wang, 1992; for deriva- 
tion). At the inner boundary, since r, > 0 and r r < If, the characteristics 
d//d< — /•,. — l , anddr/df = ty — 1/ are towards the lower boundary from interior 
(i.e., outgoing) and need to be considered. There are four incoming characteristics 
(tv, tv + l\, tv + V f , and one that is degenerate because of the model symmetries), 
so four variables can be specified at the boundary. Two other variables need to be 
calculated from related compatibility equations. We choose the values of B r , Bg, 
P' and T to be specified, leaving two quantities (i.e., tv and v 6 ) to be computed 
according to following compatibility equations: 

Oi’r _ V,B _ + V/C_ 

ot - t>v,v f (v;- - vf) ' (AJ) 

Ovg _ ^ s(^a - - V/(Vf - V£)C- 

Ot ~ V s Vf{Vj - V})B T Bg ■ (A ' 2) 

with the corresponding variables simplified in two dimensions as follows: 
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,.2 _ (# + B l) 

() / A 


(A. 2) 


V ~~ + [(tr A- lr) 2 — 4fr6“] ,/2 


70 


\ ii vvang ct a 


A TV 



r 2 


-,r _ /, 2 + [i,r + Ir 'r - 4<r/> : ] 1 - 


Z?_ - />(V/ - V;)V f (r r - l))~ - B r Bi)Vr( v r - r 

a/* 


'(‘ 7 - 1 / 


a) ('■' - - DvVftr, - r,.— 

Jr ; ,>,• 


- , Vffi 


*/■('■/ -r/nvi 

' <J 6 r </f; 

_ tWd/ - i; x -) + z?.tj -r tf r zuvr„;- — _ 

' V (W 

- L(T/ - l A 2 ) y HJ 3,,Vf ] JG, 

*• f> J /• J# /■ ,yyy 


'Wy r 


r - 2/;< r / ~ 1 a ) - ~ ~ ~!l ctg <y ) . 

flgj ctg 0 - , jf± D »BrV L _ 


1 0 v r 

~ln[Vf - r 2 , 


,>VfV * tv j- -V') 

r ■' ^ ' 


*3 77 - 

1 /■ 


• /"/0 > 2 - r 2 ; i) . 


c - - '’"'A 2 - '•?»'.«> - K)~ - B r B„Y,u-, - 


-(‘/ - >v) (iv - V,)^ + c»r, ! (iv - i: 
+ P (vl - v; 2 )v;ty-— _ 2<1Ei 

r 06 ,■ ,w 


, <>B n 

Or 


/ f uw 

~[pa z (Vl - V;) - B\) + B r B 0 V s i; 0 J- — - 

/■ Oft 

V 0 (vj- - v s 2 ) + — 4. ^ ,,r ' () b 

P J r 06 _1 " r ~of 

"y ^ 

+ ^(2*>(K 2 - V' A 2 ) 4- B 2 ] - - l Doro , 

+ ^l 


■ + 


Vift'l) , „ 

“ ( Br + Z?£ Ctg 0)4- 

“KK 2 - ^A 2 ) - B 2 j ctg # + 

-£YiYl ( V 2 _ V 2 ) _ YlEl , 


1 A. 7) 


(A. 8) 


(A. 9) 


Since the ideal MHD equations have been used, flow is parallel to the magnetic 
held lines. Thus we determine B„ from the relation B r r„ = v r B 0 . 
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